An analytic solution of the Boltzmann equation in the presence of self-generated 

magnetic fields in astrophysical plasmas 
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Through relating a self-generated magnetic field to the regular motion of free electrons that is 
responsible for the magnetic field generation in astrophysical plasmas, we solve the Boltzmann kinetic 
equation in the presence of the self-generated magnetic fields to obtain a steady-state, collisional 
invariant analytic solution of the equation. 



PACS numbers: 52.25.Dg, 95.30.Qd 

Most astrophysical plasmas are weakly inhomogeneous 
nonideal gases consisting of a large number of free 
charged particles that are closely coupled together by 
electrostatic and electromagnetic interactions. These 
interactions can propagate in the form of electrostatic 
waves such as high-frequency Langmuir waves, low- 
frequency ion-sound waves, and electromagnetic waves 
such as high- and low-frequency transverse plasma waves 
0. The so-called self-generated magnetic fields are the 
low-frequency (zero-frequency) transverse plasma waves. 
The zero frequency waves are defined to be mirror waves. 

Spontaneous excitation of the mirror waves is a well- 
established experiment fact ||] in laser-produced plasmas 
and microwave-plasma systems. In these experiments the 
free energy sources are coherent radiations. The under- 
lying mechanisms have been well understood Q. In ad- 
dition to coherent radiations, recently it has been shown 
that incoherent radiations such as thermal radiations Q 
and neutrinos |5j can also excite magnetic fields in as- 
trophysical plasmas. The keys are the nonlinear interac- 
tions bwteen plasmas and photons or neutrinos. Bing- 
ham, Dawson, Shukla, Stenflo, Bethe etc || have shown 
that nonlinear coupling between plasmas and neutrinos 
by means of the weak interaction is likely to be respon- 
sible for the most violent astrophysical activities such as 
supernova explosions. In order to take into account the 
self-generated magnetic fields in the standard stellar evo- 
lution codes we first have to solve the Boltzmann kinetic 
equation for the free electrons in the presence of the self- 
generated magnetic fields, as we are about to do in this 
letter. 

Suppose there be a self-generated magnetic field B s in 
a weakly inhomogeneous plasma. Our aim is to analyse 
how such a magnetic field affects electron distribution in 
six-dimensional phase space. For the sake of simplicity, 
we assume that the plasma consist of electrons and a 
single species of ions, the ion mass TOj be much larger 
than the electron mass m e . Real plasmas are collisional, 
so we should start from the Boltzmann kinetic equation 
of electrons ifflrol 
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where F e = (g e /m e )(E tot + v e x B tot ) is the electro- 
magnetic force, noticing that B to t contains B s . If we 
exchange subscripts e and i in Eq. (Q), we will obtain the 
Boltzmann equation for ions. 



— j\ = dp l dp' e dp' t W(p e ,p l ;p / e ,p' i ) 

■5 (i) {Pe+P l ~p' e -p' i ) 
•[/e(p;)/i(Pi)-/e(Pe)/<(Pi)]> (2) 



is the collision integral between an electron and an 
ion, in which W is the collision probability, p is a (4- 
dimensional) energy-momentum vector, while the S func- 
tion stands for the energy-momentum conservation law. 
If the subscript e is replaced by i in Eq. (||), we will obtain 
the collision integral between two electrons. 

In order to solve the kinetic equations we utilize the 
well-established two-timescale technique Q . We first de- 
compose A e = {f e , p e , E t ot,B t ot} into the fast vary- 
ing component Af and the slowly varying component 
A s , noticing that the variation time scale is similar to 
r e = 2ir/tu pe for Af, Ti — 27r/a;p C (i) for A s , where 
w pc(i) = (4'?r?T-c(i)e 2 /?7ic(i)) 1 ^ 2 is the electron (ion) plasma 
frequency. Since n 3> r e , fT' Af(t)dt = 0. Using such a 
decomposion in the Boltzmann equation and taking the 
time average over Tj , we will then obtain 
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(3) 



where F s = (? e /me)v c xB s in that E s = is assumed be- 
cause of quasineutrality in the slowly varying time scale. 
We have neglected those terms that contain ft since the 
condition ff <C f s holds good. Ion distribution function 
has only slowly varying component since an ion is much 
heavier than an electron. So we have 
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where F^ = (qi/irii)vi x B s . f s and fi couple together 
through collision integrals. 
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Generally, particle microscopic velocity v can be de- 
composed into a regular motion component u and a ran- 
dom (thermal) motion component w |?|||, 



u 



(5) 



In order to emphasize randomness of w and distinguish it 
from the regular vector u, we write w explicitly in terms 
of a random phase angle <f> and a usual vector w as follows 



i{w exp(— i(f>) + c.c.}, 



(6) 



where c.c. represents the complex conjugate of the first 
term. Obviously, w 2 = w 2 . If B s = 0, the regular mo- 
tion velocity components u e and Uj of electrons and ions 
associated with B s vanish. 

When B s is absent, F s = Fj = and we thus can apply 
Boltzmann's H-theorem M. According to this thoerem, 
we know that the steady-state solutions of Eqs. (J3J) and 
(i satisfy 



/ s (r,pi)/ s (r,p 2 ) = / s (r,pi)/ s (r,p' 2 ), 
/i(r,Pi)/i(r,p 2 ) = /i(r,pi)/i(r,p' 2 ), 
/ s (r,p e )/i(r,pi) = fs(r,p' e )fi(r,p'i), 



(7) 



which are the so-called condition of detailed balance. The 
first two have simplest general solutions as such 

ln/ s (r, p e ) = a e (r) + /3 e (r)im e [w 2 + u 2 (r)], 

In fi (r, ft) = a t (r) + ft (r) |m< [w 2 + u 2 (r)] , (8) 

since they are simplest general scalar functions made of w 
and u. The cross term u • w is meaningless and should be 
absent because of the randomness of w. The combination 
coefficient functions a and (3 depend upon macroscopic 
thermodynamic variables of the system, particle number 
density n(r) and temperature T(r), defined by 



n e (r) = -y J / s (r,p e )dp e 



T e (r) 
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and the same definitions for ions, where k B is the Boltz- 
mann constant, N e the total electron number, V the 
total volume of the system. Substituting Eq. (||) into 
these definitions, using the equality dp = mdw and the 
integration formula J Q x 2n exp(— x 2 )dx = (n — h)(n — 
I) • • • T;hy/ir, we obtain 



m e u 2 , ,n e V 



2k B l e N e 



Pe = 



ksTe 



(10) 



and the same expressions for ions. Substituting Eq. (|^) 
with ( pp| ) into the last equation of Eq. (j^), we find that 
it requires 



T e =Ti=T, 



(11) 



but imposes no constraint on regular motion velocity 
components u e and u^. 

If a single-component system is in global equilibrium, 
its particle number denisty n — N/V, regular motion 
velocity u and temperature T do not depend upon r. 
Consequently, we can make such a Lorentz boost that 
u = and a = —^ln(2irmk B T). In this special case, we 
reproduce the well-known Maxwell distribution 



/(p) - (27rmfc B T)- 3 / 2 exp 



(12) 



where w = p/m is the random thermal motion velocity 
with respect to the comoving reference system with the 
system. In a two-component system, if u e ^ u^, even if 
the system is in global equilibrium, we cannot find out 
a unique comoving reference system so that it comoves 
with both components simultaneously. Of cause, if the 
system is in only local equilibrium, there is no a unique 
comoving reference system, either. Therefore, in local 
equilibrium described by Eqs. (||), ( |Ic| ) and (pi] ), if we 
retain a = — | \n(2nmkB T), the first of Eq. (|10|) reduces 
to 
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and the same expression for ions, which are the natural 
extension from global equilibrium to local equilibrium. 
When B s is present, we use 



/ s (r,p) = (27rm e /c B T)- 3 / 2 exp 



m e (u 2 + w 2 ) 



/i(r,p) = (2n mi k B T)- z ' 2 
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as the try solution of Eqs. (||) and (Q). Substituting 
Eq. (O) into Eqs. (^) and (|j), using the explicit expres- 
sions of collision integrals, Eq. (|^), we find all collision 
integrals vanish due to the condition of detailed balance 
(0) , the terms like F • 9//<9v vanish since (wxB s )-w = 
and (u x B s ) • w — (randomness), and df jdt = 0. Con- 
sequently, Eq. (0) reduces to 



V\n(Ef es /k B T) = 0, 



(15) 



where E^ cg = ^m e u 2 is the regular motion energy per 
electron. Solving it for Ef* cg , we obtain 



K Rog = 9k R T, 



(16) 



where 9 is a dimensionless integration constant. It does 
not explicitly depend upon spatial coordinates but de- 
pends upon the magnetic field B s , the plasma tempera- 
ture T, and the number density of free electrons n e . 

In order to relate 9 to those quantities, we need ex- 
pressing Ef- eB in terms of B s first. In fact, it has been 
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shown H that the self-generated magnetic field is due to 
a slowly- varying solenoidal current given by 



167TTO e O;g 



V x (E x E*), 



(17) i£< 



where E is the slowly varying complex envelope of the 
total high-frequency electric field E: 



E = i{E(r,t)< 



— iujQt 



E*(r,i)e IW °*} 



(18) 



Inasmuch as it is the resonance between the plasma and 
radiation fields that is responsible for the magnetic field 



excitation, ujq 



the high-frequency waves can be 



considered to be quasi-monochromatic waves. So we set 
luq = uj pe hereafter. The validity condition || for Eq. Jl7j ) 
is 

k/k Dc « 1 and (|E| 2 /87rn e T e )(fc/fe De ) « 1, (19) 

where /cdc = (47m e e 2 /T e ) 1 / 2 is the electronic Debye wave 
number, k is the characteristic wave number. 

The Maxwell equation for the slowly varying electro- 
magnetic field reduces to V x B s = 47rj s in the steady- 
state limit. From this equation and Eq.([17|) we obtain 



-«(e/4m e w pe )E x E*, 



(20) 



which shows that the magnetic field self-generation in 
the plasma is a kind of nonlinear behaviours. If Eq. Jl9| ) 
is satisfied, we can neglect the convective term in the 
momentum balance equation for an electron-fluid ele- 
ment. Using this approximation and Eq. (^), we re- 
late the total fast varying electric field E to the fast 
varying electron-fluid element velocity u e , u' e (r, t) = 
— i(e/m e wo)E(r, t), where u' e is the complex envelope of 
u e , as E is the complex envelope of E. Using this relation 
and defining u' e = u r + iiij, Eq. ((2^) can be rewritten as 
follows 



B s = -(m e w pe /2e)u r x Ui 



(21) 



If we measure B s over a much larger time scale r than 
its correlation time scale tj,, B s can be considered to be 
stochastic. The mean magnetic pressure Pb over r 3> ti, 
can be expressed by 

P B =< \B s \ 2 /8n > T = {m 2 e Lo 2 pe /64ire 2 )\u r \ 2 \u % \ 2 . (22) 

Assuming u r ~ Ui without loss of generality, this equa- 
tion leads to 

£ e Rog = \m e < |u e | 2 > r = (4P s m e /n e ) 1 / 2 , (23) 

which expresses Ef cs in terms of Pb or B s . Comparing 
Eq. M) with Eq. @, we obtain 



8 = V2 



u pe ksT ' 



(24) 



which expresses the integration constant 8 in terms of the 
self-generated magnetic field B s , the plasma temperature 
T and the number density of free electrons n e , where 
Q ce = eB s /m e is the electronic cyclotron frequency with 
B s = (SttPb) 1 ^ 2 ■ 8 can be called as magnetic parameter 
since it is determined by the magnetic field B s . From 
Equations (Ell) and (M) it can be seen that 



< 8 < k De /k, 



(25) 



where = 27r/fc stands for the characteristic corre- 
lation length of the self-generated magnetic field B s . 
If Lfc is much larger than the electronic Debye length 
tdc = l/^Do as it should be, 8 can be much larger than 
unity. 

For example, using our simulation result Q B s = 6 x 
10 8 Gauss at the solar center with n e — 6.26 x 10 25 cm" 3 , 
T c = 1.36 KeV, we can estimate 6 C by means of Eq. j2^ ) 
to be 8 C = 12.5. Such large a 8 does not ruin the theory. 
In fact, numerical simulations Q show > 10~ 4 cm. 
Using T c and n e we can estimate r Do = 3x 1CU 9 cm. We 
thus know k/k Uc = 2Trr Bc /L b < 2 x 1CT 4 . Using Eq. @ 
and 8 C we can estimate (|E| 2 /87rn e T e )(/c//c Dc ) < 3x 10~ 3 , 
which shows that the validity condition ( |l9| ) of the theory 
holds good. 

Using the magnetic parameter 8, the electron distribu- 
tion in the presence of the self-generated magnetic field 
can be expressed as follows 



/ e (r,p) = / eO exp(-0), 



(26) 



where / e0 = (2Trm e kBT) 3 ^ 2 exp(— m e w 2 /2ksT) is the lo- 
cal maxwellian in that w is the random thermal motion 
velocity. 8ksT takes a role of chemical potentials, called 
as quasi-chemical potential Eq. (|2^) with 8 defined 
by Eq. (^J) is a steady-state, collisional invariant ana- 
lytic solution of the Boltamann equation for free elec- 
trons in the presence of self-generated magnetic fields in 
astrophysical plasmas that are in local thermodynamic 
equilibrium. This is the main result of the paper. 

With Eq. (120) , the statistically- averaged 7 Be electron 
capture rate |9(] in stellar interior is equal to 



A e7 = A° 7 exp(- 



(27) 



see (3.18) of [f(J for X° e7 . Using 8 C = 12.5 estimated 
above, we find the suprresion factor of the 7 Be solar neu- 
trino flux is 3.6 x 10~ 6 in the solar core since the mag- 
netic field-related regular motion of electrons will hinder 
the electron capture reaction of 7 Be. Such large a sup- 
pression factor is large enough to explain the 7 Be solar 
neutrino missing problem, the root of the solar neutrino 
problem [ pT[ . This, however, tends to worsen the 8 B so- 
lar neutrino deficit, which shows that even if there be 
self-generated magnetic fields in the solar core we still 
need matter-enhanced neutrino oscillations to fit [jllj so- 
lar neutrino experiments. Nevertheless, in this case, the 
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neutrino mixing parameters such as the neutrino mass 
and mixing angle may change. This may reduce the dif- 
ference of the neutrino masses (mixing angles) inferred 
from atmospheric neutrino experiments and solar neu- 
trino experiments. In order to investigate this important 
problem quantitatively, we have to use Eq. ( p6| ) to re- 
calculate the equation of state and opacities of the stellar 
plasma with the self-generated magnetic fields taken into 
account. We'll present such research results separately. 
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